Spin-orbit interaction in final state as possible reason 
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Abstract 



A model for ternary fission is discussed in wtiich a ttiird particle (a-particle) 
C5 I is emitted due to non-adiabatic change of the nuclear potential at neck rapture. 

An expression for energy and angular distribution of a-particles is proposed. It is 
shown that an interaction between spin of fissioning system and orbital momentum 
of a-particle (spin-orbit interaction in the final state) results in recently observed 
asymmetry of a-particle emission, which can be formally related to T-odd correla- 
^ ■ tion. No strong dependence of the asymmetry on the angle of a-particle emission 

'. with respect to the fission axis is predicted by the model in accordance with the 

experimental data. 

> ' 

^ 1 Introduction 

^ ! Recently T-odd correlation 

Cj ; B = {s [pLF X ptp]) (1) 

^ ■ was observed in ternary fission of ^ss-q g^,^^ 235u nuclei by cold polarized neutrons, where 
s is the unit vector along the neutron spin s, and plf and pxp are the unit vectors along 
the momenta plf and pxp of light fragment and a-particle (third particle), respectively 

>< : [i]-[3]. 

; One used longitudinally polarized neutrons. Let us assume that the axis y is aligned 

with the neutron beam, thus the vector s is directed with or against the axis y. The 
target is placed at the origin of the coordinates. Fragment counters are aligned so that 
they fix fragments moving with and against the axis z. Plates with a-particle counters 
are disposed bilaterally along the target and transversely to the 

Let us assume that the frame of axes is right, the vector s is aligned with the axis y, 
and the vector plf is aligned with the axis z. Therefore, when a-particle is emitted along 
the axis x, then B = +1, and B = —1 for the opposite direction of a-particle emission. 
Denote the count rate of a-particles for fixed B by N{B). Thus, the asymmetry 

iV(+l) + iV(-l) ^ ^ 

is measured. To exclude systematic errors one performs the measurements for different 
directions of the vectors s and Plf with respect to the axes y and z. The results of 
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measurements are: D = -(2.52 ± 0.14) • 10"^ for ^^^u and D = +(0.83 ± 0.11) • 10"=^ for 

235 u 

The correlation (1) is T-odd, i.e. it changes the sign along with the vectors s, p^^ 
and ptp at time reversion. However, its observation cannot be considered as an evidence 
for violation of time-reversal invariance. Similar T-odd correlation 

(S [p„ X f>LF]) , (3) 

where p„ is the unit vector along the momentum of an incident neutron, is known for a 
long time in fission of nuclei by polarized slow neutrons (see, e.g. [4]). It manifests itself 
as left-right asymmetry of light fragment emission with respect to the plane, formed by 
the vectors s and p„. 

The same left-right asymmetry is observed in elastic scattering of any transversely 
polarized particle by nuclei. This asymmetry results from the correlation 

(s[pxp']), (4) 

where s is the unit vector along the spin of the incident particle, and p and p' are the 
unit vectors along the momenta of the incident and scattered particles, respectively. The 
reason for this correlation is a spin-orbit nuclear interaction ~ si, where 1 is the relative 
orbital momentum of the incident particle and target nucleus. 

Let us assume that T-odd correlation (1) in ternary fission also results from spin-orbit 
interaction in the exit channel. This work is devoted to the theoretical analysis of this 
hypothesis. Note, that a-particle emitted in ternary fission has zero spin. Thus, we 
mean the interaction ~ Jl of spin J of the fissioning nucleus after a-particle emission 
and relative orbital momentum 1 of the a-particle and nucleus. Evidences for a spin-orbit 
interaction of such type between target nucleus spin and relative orbital momentum were 
presented, e.g. in [5]. 



2 Naive estimation 

Let us start from a naive estimation of possible effect using analogies with classical me- 
chanics and electrodynamics. As the result of a polarized neutron capture a fissioning 
nucleus is also polarized. Thus, the magnetic moment fi of the nucleus is aligned with 
the axis y. Suppose, it generates a magnetic field of the dipole type outside the nucleus 

^^Mt^^ (5) 

where n = r/r is the unit vector along the radius- vector. Therefore, the magnetic field 
is normal to the plane {x,z), in which a-particles move in the experiment [1, 2]. Let us 
assume for definiteness that the vector fi is directed against the axis y, then everywhere 
in the plane {x, z) outside the nucleus the vector H is directed along the axis y. 
The a-particle is affected by the Lorentz force 

F=-[vxH], (6) 

c 

resulting its sideways deviation (here e is the elementary charge, and v is the velocity of 
the CK-particle) . When the ct-particle moves along the axis x, it deviates in a direction of 
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the axis z, while the a-particle moving in the opposite direction deviates against the axis 
z. 

Consider the situation when a hght fragment is emitted along the axis z. It is known 
that the maximum of angular distribution of a-particle with respect to the momentum 
-pLF of the light fragment falls at the angle 9 — 82° due to stronger Coulomb repulsion 
from the heavy fragment. Therefore, in this classical picture the magnetic field deviates 
the maximum of the angular distribution of the a-particles emitted along the axis x 
toward the edge of the plate, formed by a-particle counters. At the same time for the 
a-particles emitted against the axis x the shift of the maximum of the angular distribution 
is directed toward the center of the plate of a-particle counters. Thus, one can expect 
a difference in counting of a-particles, emitted along and against the axis x, i.e. the 
observing asymmetry. 

Let us give the quantitative estimate of the effect. The magnetic field goes down as 
1/r'^, therefore, it is naturally to assume that the effect is related to a small region of 
size R, where R is the radius of the fissioning nucleus. We assume that in this region the 
velocity v of a-particle is constant. During the time At = R/v the a-particle is affected 
by a moment of the Lorentz force 

2eRH 2ea 

— - si- 

which is equal to an increase of the transversal momentum Ap_L = ma^v\_ of the a- 
particle (here = 4m is the mass of the a-particle, where m is the nucleon mass). It 
gives the angle of transversal deviation 

M^^^ (8) 
V 2mcvR^ 

The asymmetry D is of the scale of the dimensionless parameter A^. Taking for n the 
nuclear magneton e^/2mc, we get 



(9) 



he V \ mcR J 

The appearance of the fine structure constant /he = 1/137 is natural, because we con- 
sider here the electromagnetic binding between the charged a-particle and the magnetic 
field of the fissioning nucleus. 

In quantum mechanics this interaction is determined by the Hamiltonian 

(p-eA/cY 

H^^ (10) 

where A is the operator of vector potential of the electromagnetic field. In the case being 
considered it is generated by the magnetic dipole and has the form 

A = l^. (11) 

where jj, is the operator of nuclear magnetic moment. Neglecting the term ~ we 
rewrite the Hamiltonian as follows 

H=^-^t^, (12) 
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where 1 = [r x p] is the operator of orbital momentum. Taking into account that the 
magnetic moment of the nucleus is aligned with its spin, i.e. (i = fiJ/J, we find that 
the electromagnetic spin-orbit interaction ~ Jl is responsible for the sideways deviation 
of the a-particle in the magnetic field. 

But in the region of size R the nuclear spin-orbit interaction ~ Jl is at least two order 
of magnitude larger, because it does not include electromagnetic constant e^/kc. Thus, 
one can expect that the effect resulting from the nuclear spin-orbit interaction is of the 
scale 

(13) 

Taking as rough estimate R ~ 10^ h/mc ~ 10~^^ cm and v ~ 10^^ c we get the required 
magnitude D ~ 10~^. 

In a literal sense the picture proposed can be tested by studying the shifts of the 
maxima of the angular distributions of a-particles emitted along and against the axis 
X for given vectors s and Plf- These shifts would lead to a strong dependence of the 
counting rate asymmetry on the angle 9 between the momenta pxp and Plf- It was 
analyzed with the use of pairs of a-particle counters corresponding to different angles 9, 
however, no strong dependence was found (see [1, 2]). A general form for the probability 
dw of ct-particle emission into the solid angle dfl was proposed 

^r^{l + DB)F{cos9), (14) 

being in accordance with experimental data. Indeed, let us fix the vectors s and Plf 
along the axes y and z, respectively. Thus, the probability of a-particle emission in the 
direction fixed by polar 9 and azimuthal if angles is given by 

^ {1 + Dsm9cosip)F{cos9). (15) 

Here (f is the angle between the projection of the vector pxp on the {x, y) plane and the 

axis X. In the experiment one fixes the a-particles with 9 ~ 7r/2, thus taking into account 
that iV(+l) ~ dw{(f ~ 0)/dfl and A'"(— 1) ~ dw{ip ~ Tr)/dfl we obtain for the asymmetry 
(2) the quantity D. 

The classical and naive picture collapses when tested by experiment. However, it 
hardly means the failure of spin-orbit mechanism to produce T-odd correlation (1) in 
ternary fission. Note, that there is no classical limit for spin of quantum particle, therefore, 
"magnetic field", related to spin of the fissioning nucleus, cannot be treated consistently 
as classical field. This is especially true for the nuclear field of spin-orbit forces which is 
certainly non-classical. 

Thus, the problem is to construct the quantum model of ternary fission and to analyze 
in its frame a possible role of spin-orbit interaction in formation of T-odd correlation. 



3 Model for ternary fission 

The phenomenon of ct-decay is a manifestation of quantum tunneling of a-particle through 
the barrier formed by attractive nuclear and repulsive Coulomb forces. Thus, one supposes 
that the wave function of the parent nucleus includes with significant amplitude the 
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component, corresponding to two interacting by nuclear and Coulomb forces clusters - 
daughter nucleus — 4) and a-particle. 

Evidently, it is the case for fissioning nucleus also. Before scission to two fragments 
^^Al and ^"Ah the wave function of the parent nucleus ^A {Zl + Zh = Z, Al + Ah = A) 
is mainly represented by two-cluster {^^Al + ^"Ah) component. However, there is no 
reason for supression of a three-cluster component ^^A^ -\- ^"Ah + a {Zl-\- Zh — Z, 
Al + Ah + 4 = A) in the wave function of the parent nucleus ^A. 

As a result of a transience of fission process this three-cluster component cannot man- 
ifest itself via the tunneling through the barrier. But sharp changing of nuclear form and, 
consequently, of the nuclear potential can initiate non-adiabatic transition of a-particle to 
continuum from bound or quasibound state. Let us assume that such a transition results 
in ternary fission. 

The most sharp change of nuclear potential evidently takes place at neck rapture. 
Therefore, non-adiabatic transition populates continuum states described by wave func- 
tions enhanced near the neck. It corresponds to usual belief that in ternary fission a- 
particle is emitted from the neck, supported by classical trajectory calculations. 

Thus, we assume that a-particle is emitted during the short time r due to non- 
adiabatic change of nuclear potential at scission. At the final moment t = t we have the 
a-particle and two fragments ^Mx, and ^"Ah with a distance Rf between their center-of- 
masses. The origin of coordinates is in the center-of-mass of two fragments. The centers 
of fragments are located at the axis z, and this axis is directed from the heavy fragment 
to the light one. Potential energy of the ct-particle is given by the sum of potentials from 
both fragments 

Uf{r) = Ul{v - Rfi) + Uh{v + i?/z), (16) 

where z is the unit vector along the axis z. 

Really further acceleration of a-particle occurs in the time dependent Coulomb field 
of two fragments removing one from the other. To simplify the model we neglect the 
fragment movement. Then, the wave function ■0^~^(r) of the a-particle in the final state 
corresponding to the energy Ej = Ti'^k'^ /2ma is given by a solution of the stationary 
Shroedinger equation with the potential (16) with the asymptote 

p—ikr 

^^"^W /(k)-^ + e'"^ (17) 

where k is the wave vector of the a-particle at infinity. Here and below we assume for 
simplicity that the fragment Coulomb fields are screened at large distances, i.e. go down 
faster than 1/r. 

Then let us assume that at the initial moment t — Q the potential energy of a-particle 
is given by the sum of (16) and the potential Uir) caused by the neck between the 
fragments. We denote '0«(i') ^ solution of the stationary Shroedinger equation with the 
potential 

U.,{r) = Uf{r) + U{v) (18) 

for the a-particle with the energy E^. 

During the time r the Hamiltonian changes 
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Let us introduce the time dependent Hamiltonian 



H{t)=Hf + V{t). 



(20) 



The time dependent part has the hmiting values 



V{0) = U{r), 



V{t) = 0, 



(21) 



and we consider it as perturbation. Then in the first order of perturbation theory the 
probabihty of non-adiabatic transition from the initial state ipi{r) to the final continuum 
state V'k(i') is given by (see [6]) 



where Wfi = (Ef — Ei)/h. 

In the experiment one fixes the ^-particles emitted into the solid angle dVl and falling 
to the energy interval dEj. Taking the normalization volume equal unity we sum (22) 
over the final states, corresponding to the intervals dO, and dE^. Thus, we obtain for the 
differential probability of ct-particle emission into continuum 



This formula determines angular and energy distribution of a-particles in ternary fission in 
the given model. Integrating (23) over all solid angles and energies we get the probability 
of ternary fission with respect to the binary one. 

4 Scheme of calculations 

To calculate the total and differential probabilities of ternary fission given by (23) we need 
the initial and final wave functions of a-particle as well as the law V{t) determined by 
the dynamics of the final stage of fission process. 

The potentials ?7j(r) and Uf{v) acting to a-particle in the initial and final states have 
azimuthal but not spherical symmetry. Thus the projection m of the orbital momentum I 
of a-particle on the axis z conserves, but not the orbital momentum itself. 

An eigenfunction ^i(r) of the initial Hamiltonian describing a bound or quasibound 
state corresponds to the definite energy E^ and projection mj. To be more precise, there 
are two degenerated states corresponding the energy Ei and projections rrii and — m^, 
respectively. We take "01 (i") ^ ^ superposition 




T 



2 



(22) 




(23) 



V'i(r) =Vm,(r) + e'^V-m,(r), 



(24) 



normalized to unity 




(25) 
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Here 77 is a random phase. The function ipmiiT^) can be represented as the superposition 
over orbital momenta li > rrii 

^mi(r) = {Akmi) Am,(r), (r|/»m,) = i''F^,„^,(r•). (26) 
k 

Putting this function to the Shroedinger equation 

Hiil^i{v)^Eixl^i{v), (27) 
we obtain the coupled equations for radial functions 

%^ + - ^A™. + ^/..». - ^ ^(kmmi.m,)f,.^ = 0. (28) 

Together with the boundary conditions 

fkmXr) -t (29) 



1 — >oo 



they determine the initial wave function. 

Then we are looking for the wave function of the final state as a series in spherical 
harmonics 

4"^(r) = E(r|M^l;;Hr)- (30) 

Im 

Substituting it into the Shroedinger equation 

Hf^l:\:\v) = Efi:\^\v) (31) 

we get the coupled equations for the functions R\^{r) of the same type as (28). Note that 
matrix elements of potential entering the coupled equations are diagonal on the projection 
m due to azimuthal symmetry. 

To satisfy the boundary condition (17) it is convenient to take the radial functions in 
the form 

RL\r)=i:^^Y*M^I^. (32) 

lo ^ 

We get for the functions -F/^(r) corresponding to the projection m of orbital momentum 
on the axis z 

- ^-^F-ir) + k'F-ir) -^-^i:{HUf{r)\l'm)F-,{r) = 0, (33) 
with boundary conditions 

Filir) ^ {{kr)hi^\kr)6u, + {kr)h[-\kr)S^{l - k)) . (34) 

The equations (33) and (34) describe a scattering of a-particle by the deformed potential 
Uf{r) provided the ingoing state is a superposition over orbital momentum while the 
outgoing one is determined by asymptotically fixed orbital momentum Iq. 
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Substituting (26), (30) and (32) into (23) and averaging over random phase r) we get 
the double differential probability of a-particle emission as series in Legandre polynomials 

"^"^ ^ {2Q + l)aQ{Ef)PQ{cose). (35) 



dQdEf Q=o,i,2... 

The energy dependent coefficients aq^Ef) are determined by the formula 



"^^^^^ ^ i^^n^E^^-^Y ^ /l^<°oQo<°::Qo^(^om,, Ef)A*{l',m,, Ef), (36) 



lol'o 

where 



Aikm,, Ef)^Y.J (^^i'l^.^,(f)|^|/.,^.(r)i'^ll,^,(r))e^'^/^*dt. (37) 

I h 

We take into account azimuthal symmetry of the perturbation V. 

The total probability of ternary fission with respect to the binary one is given by 
integral over solid angle and energy of isotropic term in (35) 

oo 

w = An J ao{Ef)dEf, (38) 



where 



5 T-odd correlation 

As the result of slow neutron capture by a target nucleus with spin I a compound nucleus 
arises with spin Jc, where Jc = / — 1/2 or Jc = / + 1/2. If the neutron polarization is 
equal to Pn and the target nucleus is not oriented, then the polarization of compound 
nucleus with spin Jc 7^ is 



P{Jc) = < 



1 1 

- -Pn, Jc = I — 

(40) 

27 + 3 1 

■Pn, Jc = I+-- 



_ 3(27+1) 2 

The polarization axis for compound nucleus coincides with that one for captured neutron. 
Below as elsewhere above we take the axis y along the neutron beam as the polarization 
axis. 

In a binary fission the spin Jc transforms to the sum of fragment spins Jl and Jij 
and a relative orbital momentum L between two fragments. Let us assume that in a 
ternary fission the nucleus before scission emits an a-particle with a small initial orbital 
momentum 1, therefore the residual spin J can be taken approximately equal to Jc and 

p{J) ~p{Jc). 
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Thus, there is a fissioning nucleus with spin J ahgned with the axis y. If its state with 
a projection M of spin J on the axis z is described by the function ^jm, then generally 
a pure quantum state of the nucleus is given by superposition 

*j = E«m(J)*jm. (41) 

M 

A spin state of an ensemble of such nuclei is given by density matrix, averaged over the 
ensemble, namely, 

pMM'iJ) = {aM{J)a*M'{J))- (42) 
In the case being considered the density matrix takes the form 

PMM'iJ) = (^W' + 3p(J)y^^^C',%yi,(y)j . (43) 

Now let us assume that there is a nuclear interaction depending on a relative direction 
of the spin J and the orbital momentum 1 of an a-particle. Since J = + 3h + L this 
spin-orbit interaction ~ Jl may be an effective manifestation of interactions really arising 
between the angular momenta J/, and 1, in and I or L and I. Omitting the discussion 
about the nature of the spin-orbit interaction we include it as a small correction into the 
perturbation 

V = Vo+ {Vjiiil) + (ji)Vjz) . (44) 

Here the spin-orbit term is written in symmetric, evidently Hermitian form. This term 
results in changing of the spin wave function of fissioning nucleus in the process of a- 
particle emission. Therefore the model for ternary fission described above should be 
generalize to take into account the angular momenta. 

To do this note that the differential probability of ternary fission caused by pertur- 
bation (44) is given by the equation (23), where the following replacements should be 

V'i(r) — > V'i(r)*j = V'i(r)EaM(J)*jM, 

^ (45) 

4-)(r) 4-)(r)*^M,. 

We obtain the probability of a-particle emission to the intervals dVl and dEj provided the 
nuclear system remains in the state with the projection Mf of spin J on the axis z. Since 
only a-particle is registered, summation over Mj should be performed. 
It gives in linear in Vji approximation 

"^"^ {2Q + l)aQ{Ef)PQ{cos0)+p{J)cosip Yl {2Q + l)bQ{Ef)P^{cos0), 



dQdEf Q=o,i,2... Q=i,2 



(46) 
where 

Pi(cos^)=sin^^M:^ (47) 
^ cicos6' ^ ^ 

is the reduced Legandre polynomial. The coefficient bqiEf) is of the form 



bQ[Ef) - - I — , ^ 



2ln + 1 // n 



n'h'{Ef-E,r ^Q(Q + l)tipl'o + l (48) 
X Im {A{km,, Ef) {cf^^Q.B* {l',m, 1, E^) - cl^^Q-iB* (l^i -1, ^/))) , 
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where 



B{lomiX, Ef) 




dt 



J^Yi^^,)e'^f^'dt 



(49) 







Note that the expression (46) includes the needed asymmetry ~ cos ip of a-particle emis- 
sion along and against the axis x. 

In the experiment [1, 2] one fixed a-particles with ~ 7r/2 and (/? ~ or tt. Thus, 
integrating (46) over E^, we get for the angular distribution the following general form 



where D is a constant, what is close to (15). 

The asymmetry of a-particle emission at the angles </? = and </? = tt is equal to 



It is natural to expect that both functions Fi(cos6') and F2(cos6') formed by summation 
of a lot of terms, as well as their ratio, change smoothly as functions of 9. Thus, it agrees 
quahtatively with the absence of strong dependence of the effect on the angle 6 [1, 2]. 
Therefore, this experimental result cannot be considered as unambiguous evidence against 
spin-orbital mechanism of T-odd correlation formation. 

Note that the function F2(cos 9) is given by summation of the terms including reduced 

Legandre polynomials (47) which are equal to zero at ^ = 7r/2 for even indexes Q = 2, 4 

It means that the effect is due to the odd in Q terms arising as the result of mixing of states 
with even and odd orbital momenta (or even and odd parity) by the coupled equations 
(33). This mixing is in its turn caused by the absence of symmetry of Hamiltonian or, 
to be more precise, of potential (16) under space inversion because the direction of the 
axis z is fixed from the heavy fragment to the light one. Thus, in the quantum model 
being considered T-odd correlation is directly related to the asymmetry between heavy 
and light fragments as well as in the naive classical approach described above. 

6 Conclusion 

The model for ternary fission is discussed in which a third particle (a-particle) is emitted 
due to non-adiabatic change of the nuclear potential at neck rapture. An expression 
for energy and angular distribution of a-particles is proposed. Neglecting the fragment 
motion during a-particle acceleration strongly facilitates the model and does not allow to 
expect faithful reproducing of corresponding experimental data. However, one may hope 
to describe qualitatively the specific angular anisotropy of a-particles in ternary fission. 

It is shown that inclusion to the model of an interaction between spin of fissioning 
system and orbital momentum of a-particle (spin-orbit interaction in the final state) leads 




dw 



(50) 




(51) 
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to the asymmetry of a-particle emission of the same type that the observed asymmetry 
related to T-odd correlation (1). No strong dependence of the asymmetry on the angle of 
a-particle emission with respect to the fission axis is predicted by the model in accordance 
with the experimental data. 

I am grateful for helpful discussions to V.E.Bunakov, W.I.Furman, F.Goenncnwcin, 
G.A.Petrov and G.V.Danilyan. The work is supported by RFBR grant 00-15-96590 and 
INTAS grant 99-0229. 
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